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Abstract 

The correlation functions of the Z-invariant Ising model are calculated explicitely 
using the Vertex Operators language developed by the Kyoto school. 
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Introduction In the last few years the Kyoto school has made a breakthrough in the 
calculation of correlation functions and form factors of Integrable Models on the lattice by 
treating them directly in the thermodynamic limit. This was possible because the physical 
quantities (observables and space of states) could be translated into the powerful language 
of representation theory of affine quantum algebras. Fundamental for this progress was an 
understanding of the combinatorial aspects of the calculation of order parameters through 
Baxter's corner transfer matrix formalism. The next step was the interpretation of the CTM 
Hamiltonian as the derivation generator in the quantum affine algebra. Finally the complete 
apparatus of representation theory developed in was employed with the introduction of 
intertwining operators for different representations - the so-called Vertex Operators. 

The success of this method has stimulated the search for new deformed affine-algebras 
and their representation theory. It has also been possible to calculate correlation functions 
without the knowledge of the dynamical symmetry in a way similar to the bootstrap 
program for form-factors in integrable field theory. Interesting also is the development in 
|r3| where techniques of Conformal Field Theory were directly 'deformed' to the setting of 
RSOS models. This progress proves interesting not only in understanding the integrability 
of lattice models but also around in understanding the issue of 'quantum symmetry' in CFT 
where up to know the concept of quantum group has not been of much use and where the 
more abstract methods of algebraic quantum field theory have yet to ripen. For this reason 
the study of the simplest integrable model, the Ising model, with this new method proves 
to be a good laboratory for understanding rigorously the new methods independent of the 
setting of quantum algebras. 

In this letter we restrict to the computation of physical quantities and show how the 
methods developed in are surprisingly simple in obtaining correlation functions that 
by the old transfer matrix methods are technically intricate. In the above paper it was 
shown how to calculate the nearest diagonal order-order and disorder-disorder correlations 
functions; we extend this without further pain to a general distance and show by using a 
theta functions identity the validity of an old conjecture of Baxter but we can go further 
and integrate this equation obtaining an explicit expression for the two-point correlation 
functions. In the last paragraph we make some comments on futher work possible using the 
methods of this letter. 
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Vertex Operators We restrict to the ferromagnetic Ising model in the edge description, 
below the critical temperature i.e. with < k < 1, k = (sinh2i^ sinh2L)^^. Here i^' > and 
L > are the coupling constants in the two orthogonal directions and we use Onsanger's 
elliptic parameterization: sinh(2ii') = —isn{iu), sinh(2L) = i/ksn{iu) where sn(-u) is a 
Jacobian elliptic function with half-periods /, il', k is the modulus, and < u < I' (see 
e.g. [0]). We will also employ the multiplicative parameters ( := exp(— 7rn/2/) and x : = 
exp(-7r/72/)). 

Viewing now the two-dimensional lattice in the diagonal and dividing the lattice in four 
equal quadrants to apply Baxter's CTM method we have two cases according to if the central 
point belongs or not to the lattice; this defines two CTM's that act on the two sectors of the 
model, NS and R, respectively (in analogy to CFT). Diagonalising the adjoint action of the 
CTM Hamiltonians we can describe the two sectors as Fock spaces (called Ti.^^ and 7i^) of 
two irreducible representations of free fermions. We have: 

[D^^0f] = -2r0f, [D\<pf] = -2r<pf, 

where r e Z + | for the NS sector, r G Z for the R sector and D^^, are the respective 
CTM hamiltonians. 

There are two vertex operators (VOs in the following) depending on which sector they 
intertwine: 

*^s(C):^^' — ^^ 

$NS(^)^^R^^NS^ 

accordingly they intertwine the above fermions in a manner which can be written in the 
surprising elegant way as follows: 

0^S(z)<i>-(C) = /«^)$-'^(C)0^(^), 
0^(^)<|>.(C) = /«^)$-.(C)0^'(^), (1) 

a$-(C) = <|.'^(C)^f(C), 

a$.(C) = -^$_.(CX^(C). (2) 
Here we use the simplified notation of g $r '"(C) = $"(C), 

a being equal to +1 or — 1 

distinguishes the two sub-sectors (with even or odd parity) of the NS sector and $nS(t(C) = 

2 



<l>(j(C), also we define [z] := f{z) = —^/ksn{v) with z = exp{i7TV / 1) and we have introduced 
the generating functions 



[0^(z), (f{w)]+ = 6^{x^zlw) + 



and 



where 



[2;Jns = \/ cn('i; 



rez+i 



rjr = X + X 

By crossing symmetry i.e. m — m or in the multiplicative language C — > we can 

obtain (using also the transformation f{x'^z) = l/f{z)) similar relations with the parameter 
of the VO replaced with x(. Finally the VOs satisfy the following unitary relations: 

J2^^i<WiO = 9''^^dnn, (4) 
<l>'^(xC)<f.(C)=^7'''xid«Ns.., (5) 

where the scalars and 5''"^'^ are given by 

4. ^4 ^8\2 f^8. ^4 ^8\2 



R l-*- ) J oo „NS I^X ,X ,X ; 



oo 



('^,2. ^,4 ^,4"! ' ^ /'t,6- ^,4 ^4"! 

For the infinite product notations see e.g. [ffl. 
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Two-point correlation functions We can now calculate the correlation functions (CF) 
for the model with an arbitrary number of inhomogeneities. As explained in [Q the CF are 
giving as traces of the VOs, in an Appendix of the same paper the nearest order-order and 
disorder- disorder CF were obtained. The generalization to arbitrary separation is immediate. 
Let us call them 5'2n('^i' "^2, • • • , C(2n) and g2j^{ai, 0^2, . . . , a2n) respectively (when we don't 
want to emphasize the spectral parameters and g!^^ simply), where (j = exp{i7Taj/2I), 
then we have: 

Ea-M. ^' f^f^'trNS (^x^^'''' F{a, a'] /Xi; Ci, • • • , C2n)) 

,2-„(a„ a,,..., a2„) = (g^SgK^n,,^^ (^2^NS) , (6) 

F(a,(T';/ii;Ci,...,C2n) = 

$"'K2n)$M„_,K2n-l) ' " " $'^n<2) Kl) (Cl ) (C2) " ' " $^"-KC2n-l ) (C2n) , 

where the subscript /i, in the sum above means the sum for all 1 < i < n — 1. Similarly for 
the disorder- disorder case: 

E^r (x'^^'^Gifii; Ci, • • • , C2n)) 
,2^Ja„ a,,..., a,.) = (^ns.R).^^^ , (7) 

Gifii] Cl; • • • ) C2n) = 

$M„(xC2n)$^"(xC2n-l) " " " ^2) -^--Mi (0) $ "^^(2) " " " (Csn-l) (Csn) , 

with 1 < i < n. To obtain explicit expressions for the above we need the following equalities 
that can be proved using eqs. (|l|),(|]), (H) and (||) (see also (Q): 

Y,^aMl)^'''i<^)^-'^AC^)^"''iCl) = ^9''^ / / ^/jR [w;.^]ns N^] 0^(^z)0^(t/^.) , (8) 

"1 

5^a2„<l>-'^^"(a;C2n)<fa.„(C2n) = -^^?''' / [W2n2n]^s<p''^{w2n), (9) 

where we changed the variables to Wi = ZijQ and used the notation §^ = § 2^7^; ["W^jj] = 
Because of eq. (||) the fermions in expression (H) can be anticommutated with 



only a sign change. In the case of (^gn need only relation starting at the center 
in the numerator of expression (|^) we interchange each time four VOs with to fermions, 
these we intertwine to the extreme left or right using (|1]) paying attention to the change of 
representations we do this successively n times to arrive at: 

» 2n 
l<i<2n-^^i l<i<j<2n i=l 

where we have 2n integrations and the product of fermions is ordered from left to right with 
increasing indices. Similarly for the (^gn but starting with the first eq. in @ and ending 
with the second one, the final expression is obtained interchanging the R with NS in the 
expression above. At the end we are left with, using the notation 6^{ij) = 6^{x'^Wi/wj), 

tr„B,(x2^>I^(;.O0^(^2)---0^(^2n-l)0«(;^2n)) ^, ^ 

^ : =E iPiP2)S iPm) ■ ■■^''iP^.-.Pj- 

llfiB [X ) ^ 

Here pi, . . . ,P2n is some permutation of the numbers 1, 2, ... , 2n, Yl'p is a summation over 
all {2n — 1)!! permutations which satisfy the restrictions P2m-i < P2m for 1 < m < n and 
P2m-i < P2m+i foT 1 < 171 < Ti — 1 and Ep is the parity of the permutation. There is a similar 
expression for the NS case (note that the integrations in this case are well defined although 
the function S^^{ij) has half- integral powers in the integration variable z) . 

Before inserting this formula under the integral expression it is more convenient to an- 
tisymmetrize the factor multiplying the fermions in all the Wi variables (under the integrals 
the fermions can be anticommutated as already observed), for this we change the spectral 
parameters and integration variables to their additive version, 

dvi 



^^S'^(zj)^2/5^s>,-t;,+^/'), 



-/ 2/ 



i2Ik'}li(vi-iai) , 2Ik' Qi{vi-ia 



7T &{vi - iai) ' V vr &{vi - iui) ' 

where we introduced the Jacobi theta functions H,6,Hi and 61 (see e.g. [||). The result of 
the antisymmetrization is: 

nr=i Hi(v2i-i - ia2i-i)Qi{v2i - ia2i) ni<i<j<2n H(^^i - iaj) 



+ perm 



ni<.<,<2n0K-^«.)H(t^.-^,)(Hi(o)ei(o)r . 

- xHi 2^t;i -«2^aJ , (10) 



Hi(o)n-:in-:i0(^.-^«,o 
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A=l r=l 



where by perm we mean the antisymmetrization in all Vi. The same formula is valid under 
the interchange of 61 with Hi. This identity is proved showing first that both sides transform 
identically under the periods 21 and 2il' in all Vi, and then through an inductive process it 
is possible show that the residues of the simple poles at e.g. Vi = iai + iV are identical for 
the two sides. As a simple consequence of Liouville's theorem the rhs can only differ by an 
additive constant, which is proved to be zero by comparing the two sides at e.g. = ia2n- 
Substituting all these expressions and integrating out the delta functions, we obtain: 

fl'2n.(«ir--,«2n) = 

1 /-fc2O'^(0)\"-p-i- r l<r<s<2n l<j<k<n 

n 

:vl>„(2^t;,-z^a,), (11) 



n 2n 
j=l r=l 



Hi(u) if n even 
0i(u) if n odd 



with h{u) = 0(n)H(n). This formula proves Baxter's conjecture in (the formula was 
proved in the critical case in [0]). We obtain ^fgn interchanging Bi with Hi in ([Tl|). 

We can go further and combine the above theta function identity with the following one 
(this is a shifted version of Fay's identity, see e.g. [Q): 



Hi(t;,-m,)H(0)'^ .TT,.n^^2n 11 K ^ jJ K ^ 



l<i<j<2n 



Hi(o) n0(^.-^«,) '[^^' 



and another identity interchanging Hi with Gi. We see that multiplying the rhs of the above 
eq. by appropriate constants terms we obtain the rhs of ([T0|) , expanding the determinant 
we have reduced the problem to independent known quadratures, collecting all constants we 
arrive at the final result: 



Pfafiian {h^''^{ai — OLj)) 



PfafRan f yksn{iai — iaj^ 
, 21 fk'Y Oi'iia) , 2/,,,,iHi'(za) 
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Here we also employed a lemma proved in |T^ (eq. (5.8)): 

Pfaffian (^\/ksn{iai — iaj)^ = Y\ (^y/ksn{iai — iaj] 

l<i<j<2n 

Expression (|T2|) is symmetric in all spectral parameters and satisfies the recurrence rela- 
tion of Baxter 

92ni.^l^ "2, ■ ■ ■ , a2n-l, a2n-l + I') = fi'2n-2("l' "2, " " " , "2^-2), 

+ /') = 1. 

As a final check of the expression for the correlation functions we calculated the homogeneous 
limit, when all spectral parameters are equal of the correlation function, taking the 
necessary limits we obtain after a tedious calculation: 

giHomo. = ^[(5 - k')E{kr + 8{k' - i)E{k)i{k) + 3{e - irm'] 

where I{k), E{k) are the complete elliptic integrals of the first and second kinds. This result 
coincides with the value obtained by Toeplitz determinant method (see [0). 

Employing similar methods as for the two-point functions we can write immediately the 



higher CP as a multiple integral expression (see ||T5[|). 



Conclusions We have shown that treating the Ising model at the thermodynamic limit 
simplifies considerably the calculations of physical quantities (we remark that using the full 
Z-invariance of the model gives us simpler results than the homogeneous case). It might 
also be possible to use the theta function identities to explicitly integrate the correlation 
functions of the six- vertex model obtained in [|10] (see also |]ll|) through bosonization. Using 



the commutation relations of the VO and its unitary properties we can write difference 
equations, it will be interesting to understand in this way the CP , this is important to 
understand the eight-vertex model CP out of its decoupling point (where it decomposes in 
two Ising models) where we cannot use the translation to fermions but the problem can also 
be formulated as the solutions of some difference equations as in the form-factor program of 
integrable field theory. The simplicity of the result ( |T2|) gives some hope in the possibility 
of generalizing this result for the eight-vertex model (by some deformation with an extra 
parameter). Purther work is also suggested by [O where a deformed energy tensor field 



plays a role of dynamical symmetry for the RSOS models, the Ising being the simplest one 
in this category it will be interesting to understand in this context this symmetry and its 
connection or not with the temperature-like states that happen to describe the CF of the 
theory and the toroidal parameter space of the elliptic functions used (see 0). 
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